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Abstract 

Using the Dirac-Brueckner-Hartree-Fock approach, the properties of neutron-star matter includ¬ 
ing hyperons are investigated. In the calculation, we consider both time and space components of 
the vector self-energies of baryons as well as the scalar ones. Furthermore, the effect of negative- 
energy states of baryons is partly taken into account. We obtain the maximum neutron-star mass 
of 2.08 Mq, which is consistent with the recently observed, massive neutron stars. We discuss a 
universal, repulsive three-body force for hyperons in matter. 
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Neutron stars may be the most dense and exotic state of nuclear matter, and its core 
serves as a natural laboratory to investigate the nuclear matter whose density reaches several 
times higher than the normal nuclear-matter density, n° B [1]. In fact, the recently observed, 
massive neutron stars, J1614-2230 (the mass of 1.97 ±0.04M e , M 0 : the solar mass) [2] and 
J0348+0432 (2.01 ±O.O4M 0 ) [3], have provided important information on the equation of 
state (EoS) for dense nuclear matter. 

To understand these heavy objects, various nuclear models have been examined, in which 
relativistic mean-field theory (RMFT) is very popular and has been successfully applied 
to the dense nuclear matter [4], However, in RMFT, nucleon (N)-nucleon short-range cor¬ 
relations in matter cannot be treated. In contrast, in the Dirac-Brueckner-Hartree-Fock 
(DBHF) approach, although the calculation is involved, one can consider the effects of the 
Pauli exclusion principle and short-range correlations. 

Until now, several groups have performed the DBHF calculations not only in the region 
around n° B but also in matter at higher densities (see Refs. [5-15]). However, so far there 
has not been any relativistic attempt to take account of the degrees of freedom of hyperons 
(Ys) as well as nucleons in dense matter. Because it is quite interesting to see how hyperons 
contribute to the EoS and to the maximum mass of neutron stars, it seems very urgent to 
perform the DBHF calculation for matter including hyperons. 

In this Letter, we study such dense neutron-star matter using the DBHF approach. Here, 
we particularly pay attention to the following two points: (1) the space component of vector 
self-energy of baryon (B), is taken into account, because, although it is certainly small 
at low density, it is expected to be important in dense matter, (2) as in Refs. [13-15], we 
partly consider the effect of negative-energy states of baryons in the Bethe-Salpeter (BS) 
equation to remove the ambiguity in the relationship between the on-shell T-matrix for 
baryon-baryon scattering and the baryon self-energies [8-10]. Furthermore, when hyperons 
take place in matter, the effective masses of interacting two baryons become very different 
from each other, and thus we should treat the baryon-mass difference in the BS equation 
explicitly. 

We now start with the self-energy of baryon in the rest frame of infinite, uniform nuclear 
matter. It is given by 

£ B (k) = Ef (k) - To E° (*0 + 7 ■ kS v B (k), (1) 
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where k (. k ) is the three (four) momentum of baryon. Here, T, B ( ' 0 ' 1 ^ is the scalar (zero-th 
component of vector) [space component of vector] part of baryon self-energy. Using these 
self-energies, the effective mass, M B , the effective momentum, k* B , and the effective energy, 
E B , in matter are defined by 

M- B (k) = M B + T. s B (k), k' B = fc[l + Sg(fc)], E' B (k) = Jkg + Mg{k), (2) 

with Mb being the free baryon mass. Then, the baryon spinor states with positive or 
negative energy are respectively constructed as 

<mm) = ^Mm+EE) ( ), ( 3 ) 

\ M*(k)+E* B (k)Xs ) 

e B (k, s) = yj M’ B (k) + E%(k) ^ j , (4) 

where cr is the Pauli matrix, and denotes a 2-component Pauli spinor. 

In the conventional DBHF calculation, the baryon-baryon scattering is usually evaluated 
in the center of mass frame with respect to the interacting two baryons. In such cases, 
instead of Eqs.(3)-(4), the hclicity spinors and the partial-wave decomposition are often 
used to solve the BS equation [5-11, 13-15]. However, when T} b remains finite and k ^ k* B , 
although k and k* B are parallel with each other in the nuclear-matter rest frame, they are 

not in the center of mass frame. It is thus more convenient to perform the calculation with 

the standard spinors, Eqs.(3)-(4), in the nuclear-matter rest frame, rather than with the 
helicity spinors in the center of mass frame. 

Furthermore, the inclusion of negative-energy states of baryon in the BS amplitude may 
be necessary to remove the ambiguity of the relationship between the reaction matrices for 
baryon-baryon scattering and the baryon self-energies [13-15]. Thus, we here define four 
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reaction amplitudes 


Tb'"B"B'b {}^, k, s'", s", s', s ; P) 

= t ip + k\ s'd 4 b „ tip - k\ A r «> B , t ip + k, A #B tip - k, d , (5) 

Rb'" B" B' B^ , k, s'", s", s', s; P) 

= e B >» Qp + k', s"^j Qp - s") r $ B , Qp + k, /) Qp - k, s^j, (6) 

Ob'"b"B'B{ k', k, s'", s", s', s; P) 

P + k', s'"^j § B " 

PB'"B"B'B{k', k, s'", s", s', s; P ) 

= 0B- Q 

where T represents the effective reaction operator, and these amplitudes satisfy the following, 
coupled BS equations 


P + k', s'" $ 


■ B" 


-P 

2 


k',s”) re 


>B' 


i P + k,s' 


4, 


-P 

2 


Ml, (8) 



P — k', s" ) re B , -P + fc, S 'U B l-p-k,s , (7) 


Tbb'bb' {k, k, s, s', s, s'; P ) = Vbb'bb' ( k,k,s , s', s, s'; P ) 

+ l 7o~ \/^ BB'B"B'"{k, q, s, s', s", s'"; P) 

s"s'"B"B'" J ' 77 ' 

xQb"B"'(P , q)gThB"B'"{P , q)Tb"B'"bb'(q, k, s'", s", s, s'; P ), (9) 

RBB'BB'{k, k, s, s', s, s'; P) = U B B'BB'(k , k, s, s', s, s'; P ) 

+ / 7^31 U B B'B"B'"{k,q,s,s',s",s"';P) 

s"s'"B"B"' ' 

xQB"B"'(P,q)g ThB "B'"(P,q)TB"B'"BB'{q, k, s'”, s", s, s'; P ), (10) 

0 B B'BB'(k, k, s, s', s, s'; P) = W BB 'BB'(k, k, s, s', s, s'; P) 

+ [ 7?—A ^ BB 'B"B’"{k,q,s,s' ,s",s'";P) 

s"s'"B" B'" J ' 71 ' 

xQb"B"'{P ? q ) 9Th B" B'" ( Pi q)OB"B'"BB'{q ! fe) S W , s", S, s '; P), (11) 

Pbb'bb' (k, k, s, s', s, s'; P) = Z BB 'BB'(k, k, s, s', s, s'; P) 

+ f 7 ^AjUBB'B"B'"(k, q,s,s', s", s"';P) 

s"s’"B"B'" J ' 71 ’ 

xQb"B"'(P, q)9ThB"B'"(P, q)OB"B'"BB’{q 1 k, s'", s", s, s'; P), (12) 


with V, U, W and Z being the anti-symmetrized matrices of one-boson-exchange (OBE) 
interaction [16] with respect to the positive- and negative-energy states (as seen in Eqs.(5)- 
( 8 ))- 
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In Eqs.(9)-(12), Qbb' is the Pauli exclusion operator for baryons B and B', and grhBB' 
denotes the Thompson’s two-particle propagator [17]. The seven arguments in the four re¬ 
action amplitudes, T, R , O, P , are as follows: from left to right, the first variable represents 
the final (or intermediate) relative three-momentum; the second, the initial (or intermediate) 
relative three-momentum; the third and fourth are for the spins of the final (or intermediate) 
two baryons, each of which is up (+) or down (—); the fifth and sixth, the spins of the initial 
(or intermediate) two baryons; and the last one is the total three-momentum of interacting 
two baryons. We note that the negative-energy states appear only in the initial and/or final 
states of the BS amplitudes, and they are not included in the intermediate states, because, 
in the realistic baryon-baryon potentials such as the Bonn potentials, the negative-energy 
states are usually not considered [13, 14]. 

The ladder-approximated, coupled BS equations can be numerically solved in the nuclear- 
matter rest frame. To reduce the number of variables and make the present calculation 
feasible, we here average the azimuthal angle in the spinors, Eqs.(3)-(4), namely we replace 
E B (l/2P±k) by the averaged one, ^ / d(pE B (l/2P±k). We have checked that this change 
does not lead any large numerical error in our final results. 

Given the reaction amplitudes, we can calculate the following components [14] 


£f M = $ B (k,+)E B (k)$ B (k,+) 

= 2M* B (k)E s B (k ) - 2 E* B (k)E° B (k) + 2k ■ k* B Y^{k), (13) 

£f *(fc) = e B (k,+)E B (k)^> B (k,-) 

= 2\k B \(k)E° B (k) - 2\k\E* B (k)E v B (k), (14) 

£f e (k) = e B (k,+)E B (k)Q B (k,+) 

= - 2M* B (k)E s B (k ) - 2E* B (k)E° B {k ) + 2k • k* B E v B {k ), (15) 


where E B {k) is given by Eq.(l), and the components, Sf $ (/c), Eq $ (A;), Eq 0 (/c), are respec¬ 
tively calculated through the reaction amplitudes, T, R , P. Using these relations, we can 
uniquely determine the self-energies in Eq.(l), and calculate the energy density and pressure 
of matter [7]. We here discard the contribution of retardation effect. 

Now we are in a position to show our results. Through the whole calculation, we adopt 
the Bonn potentials [5], and use the conventional “reference spectrum” approximation [5- 
12, 15], where the momentum dependence of the self-energies is frozen at some reference 
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TABLE I: Calculated properties of symmetric nuclear matter at the saturation point, n B . In the 
first column, the results of Bonn A, B and C are respectively labeled by A, B and C, while A* 
denotes the result of Bonn A with the modified coupling g* NN(J (see Eq.(16)). The values of the 
binding energy per particle, S/n° B — Mjv, the incompressibility, K, the symmetry energy, S, and 
the slope parameter, L, are in MeV, and n^, is in fm~ 3 . The fitted values in A* are denoted by f. 


case 

„o 

n B 

S/n° B — Afjv 

K 

S 

L 

A 

0.149 

-10.5 

204 

28.8 

78.6 

B 

0.130 

-7.3 

133 

22.7 

58.2 

C 

0.112 

-5.2 

87 

18.0 

42.2 

A* 

0.1681 

-15.31 

2331 

33.61 

95.0 


momentum. Here, the reference point is chosen to be the Fermi momentum, kp, at each 
nuclear density, rip- 

We first study the symmetric nuclear matter around n ° B , where it consists of only nucleons 
interacting through the exchanges of cr, 5, u, p, r) and tt mesons. In Table I, we present the 
properties of matter at n° B . In the present calculation, the binding per particle in A~C 
is relatively shallower than the empirical value. This tendency is close to the result by 
Poschenrieder and Weigel [13] because our method resembles their approach. 

We try to adjust the matter properties by assuming that the nucleon-a coupling constant 
varies as a function of the scalar self-energy £]y. Because baryon is a composite object, 
the meson-baryon coupling strength may generally depend on the scalar density in matter 
[15, 18-20]. To take account of such an effect, we suppose that the coupling is expressed as 



where g^Na is the value in vacuum. The four parameters, a i= 1 ^ 4 , are determined so as to 
reproduce the empirical binding value, n B , K and S (see Table I), and we then find an = 
—0.36, a 2 = —1-69, «3 = —3.07 and 0.4 = —1.86 for Bonn A potential. This modification 
enhances g^Na by only 2% at n° B . The properties of matter in this scheme (A*) is also shown 
in the table. 

We notice that, in the present calculation, the space component of vector self-energy, £j^, 
is certainly small around n Bl but it grows rapidly at high density and reaches about 0.7 at 
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ns = 1.0 fm -3 in symmetric nuclear matter (see Eqs.(l) and (2)). Therefore, it is necessary 
to take account of the space part explicitly in dense matter. 

Next we challenge the calculation of neutron-star matter including hyperons. In the 
following calculations, we adopt the scheme A*. 

We now have to determine the coupling constants for hyperons. Using the experimen¬ 
tal data of nucleon-hyperon scattering, the hyperon-meson coupling constants have been 
studied by several groups [21]. However, due to poor experimental accuracy, the coupling 
constants cannot be determined without large ambiguities. Thus, in the present calculation, 
we determine them (except for the hyperon-cr couplings) by using SU(6) symmetry [20, 22], 
For the hyperon-cr interactions, we can use the recent analyses of hypernuclei and hyperon 
production reactions, which have suggested that the A, E~ and 5“ may respectively feel the 
optical potential, Ua-,e-,e- — —27, +30, —15 MeV, in nuclear matter [23]. We thus choose 
the coupling constants, gw a, so as to reproduce these potential depths at n ° B , using the 
Schroedinger-equivalent optical potential 

U Y (k) = E s y(k) - 5Mh[ S . (fc) - (fe)] + +-[(Ef(*:)) 2 - (E» (*)) 2 ]. (17) 

Furthermore, for nonstrange mesons, a cutoff parameter in the form factor at hyperon- 
meson vertex, A yy'm, is assumed to be the same as in the nucleon-meson form factor, while, 
for strange mesons, a cutoff parameter, A bb'k (A bb'k*), is taken to be the average value of 
A bb’t] and A bb’-k (A bb'u and A bb'p)- 

Because an enormous amount of time is needed to perform the full calculation, we first 
perform two preliminary calculations: one includes e ~, p~ and five baryons (neutron (n), pro¬ 
ton (p), A, E _ , S - ), 1 and the other includes the leptons and eight baryons (n, p, A, E _ , E°, 
E + , S - , S°). In these calculations, we consider only six nonstrange mesons (a, 5, oj, p, rj, 7r), 
and exclude the baryon-exchange and baryon-transition processes such as N + A —> A + N, 
N + A —* N + E, etc. Note that these processes are induced by the exchanges of K, K* and 
iso-vector, nonstrange mesons. We call the first set B5M6, and the second B8M6. 

In the panels (a) and (b) of Fig.l, we present the particle fractions in B5M6 and B8M6. 
As seen in the figure, both the results are very similar to each other, and show that the 

1 From among the members of the E and H hyperons, we select the E _ and only. The reason is because, 
from the viewpoint of electric charge, it is expected that they can appear easier in matter rather than the 
other members (E + , E° and S°) [20]. 
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FIG. 1: (Color online) Particle fractions for (a) B5M6, (b) B8M6 and (c) B5M8. 

X - first appears around % ~ 0.37 fm -3 , and next the S - and A are created. However, the 
fraction of XA dwindles rapidly with increasing ns- In contrast, the numbers of 5“ and A 
grow steadily, once they emerge in matter. 

Because the difference between B5M6 and B 8 M 6 is expected to be small, we proceed 
to the final calculation, where the five baryons are considered and they interact through 
the exchanges of eight mesons (a, 5, u>, p , p, 7 r, K, K*). We here include the effect of the 
baryon-exchange and baryon-transition processes. We call this scheme B5M8. In the panel 
(c) of Fig.l, the particle fraction for B5M8 is displayed. It is interesting to notice that 
the result is again similar to the previous ones, but the fraction of X - is enhanced by the 
baryon-transition process between X - and A hyperons. Furthermore, comparing with the 
result in B5M6, the thresholds for the A and S _ in B5M8 move toward higher density. 

Using the Tolman-Oppenheimer-Volkoff (TOY) equation [24] with the BPS model [25] for 











R (km) 


FIG. 2: (Color online) Mass-radius relations for neutron stars. The dot on each line represents the 
maximum mass (see also Table II). The shaded area represents the mass of J0348+0432. 


TABLE II: Neutron-star radius, R m ax (in km), the central density, n c (in fm 3 ), and the ratio of 
the maximum neutron-star mass to the solar mass, M max /M@. 


case 

Rmax 

n c 

M m ax /Mq 

NM6 

12.6 

0.78 

2.44 

B5M6 

13.1 

0.76 

2.03 

B8M6 

13.1 

0.75 

2.04 

B5M8 

13.1 

0.74 

2.08 


the EoS in the crust region, we can calculate the neutron-star mass as a function of its radius. 
The calculation is performed under the conditions of charge neutrality and /3-equilibrium 
in weak interaction. The present results are summarized in Fig.2 and Table II, where we 
show the mass-radius relations and the properties of neutron stars at the maximum mass. 
Here, NM6 denotes the result in which only the leptons, nucleons and six nonstrange mesons 
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participate. We can find that the predicted maximum mass in each case is consistent with 
the observed ones, 1.97 ± 0.04 M 0 (J1614-2230) [2] and 2.01 ± 0.04 M & (J0348+0432) [3], 

We would like to compare the present results with those in the non-relativistic Brueckner- 
Hartree-Fock (BHF) approach. It is well recognized that the saturation properties of sym¬ 
metric nuclear matter cannot be explained by the BHF calculation with two-body inter¬ 
actions, and that it is vital to consider a repulsive three-body force (TBF) additionally 
to move the calculated saturation point toward the empirical value [26]. In contrast, the 
DBHF calculation can provide a result close to the empirical value without further ingre¬ 
dients. This is because the DBHF approach involves an inherent ability to account for 
important TBFs through its density dependence, i.e., the TBF originating from virtual ex¬ 
citation of a nucleon-antinucleon pair, known as Z-graphs [5]. In fact, it is not difficult to 
find the term of Z-graphs, A E^ air rj k 2 (S^) 2 / 2M %, in the nucleon energy at the mean-field 
level [27-29]. The relativistic effect on the binding energy per nucleon is then well fitted 
as A (S/n%) re i oc (ng/n^) 8 / 3 [5], which depends on rig strongly, and it helps obtain a bet¬ 
ter saturation point [30]. From this point of view, the contribution of Z-graphs plays an 
important role in success of Dirac phenomenology. On the other hand, some people have 
argued that the pair creation should be suppressed by the compositeness of nucleon [27, 31]. 2 
However, even when the effect of the compositeness is taken into account, the repulsive term 
still remains [18, 29]. 

In dense neutron-star matter, where hyperons can also participate, the similar situation 
may occur. In this case, it is again well known that, in the BHF approach, the inclusion of 
hyperons softens the EoS very much, and that such an EoS is inconsistent with the existence 
of heavy neutron stars [32]. To remedy this problem, it may be again necessary to introduce 
repulsive TBFs for hyperons [33]. On the other hand, the DBHF calculation inherently 
contains the density-dependent, repulsive TBF, and it seems to be universal for all baryons 
[18, 29]. As long as the magnitudes of E'| and are large, due to Lorentz structure, each 
baryon feels the repulsive potential, A E^ airl in nuclear matter, which may again play an 
important part in obtaining the EoS for sustaining the massive neutron stars, as shown in 
the present calculation. In Fig.3, we show the self-energies for hyperons as well as nucleons. 

2 In Ref. [27], it has been emphasized that the “Z-graphs” in Dirac phenomenology should not be interpreted 
as arising from virtual nucleon-antinucleon pairs. 
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We can see that their magnitudes are certainly of the order of lOOMeV at high densities. 



n B ( fm ' 3 ) 


FIG. 3: (Color online) Self-energies, T, B and T,° B , in B5M8. 

In summary, using the Dirac-Brueckner-Hartree-Fock approach, we have studied the prop¬ 
erties of neutron-star matter including hyperons. The result has shown that the FT, A and 
5“ appear in dense matter, but the fraction of FA decreases with increasing ub■ The max¬ 
imum neutron-star mass is estimated to be 2.08 M e , which is consistent with the masses 
of heavy neutron stars. Thus, we can conclude that it is very important to consider not 
only the effects of Pauli exclusion principle and short-range correlations in matter but also 
the relativistic effect involved in Dirac phenomenology. In the present calculation, we have 
determined the hyperon-meson couplings by SU(6) symmetry. However, those couplings 
should be improved in the future calculation. In the DBHF approach, it is very difficult to 
understand a neutron star with heavier mass than that of J0348+0432. Thus, if such an 
object is found in the future, it may clearly suggest that the nonbaryonic degrees of freedom 
(like quarks) emerge in its core region. 
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